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Recently, several hybrid approaches to quantum information emerged which utilize both continuous- and 
discrete-variable methods and resources at the same time. In this work, we investigate the bipartite hybrid 
entanglement between a finite-dimensional, discrete-variable quantum system and an infinite-dimensional, 
continuous-variable quantum system. A classification scheme is presented leading to a distinction between 
pure hybrid entangled states, mixed hybrid entangled states (those effectively supported by an overall finite- 
dimensional Hilbert space), and so-called truly hybrid entangled states (those which cannot be described in 
an overall finite-dimensional Hilbert space). Examples for states of each regime are given and entanglement 
witnessing as well as quantification are discussed. In particular, using the channel map of a thermal photon 
noise channel, we find that true hybrid entanglement naturally occurs in physically important settings. Finally, 
extensions from bipartite to multipartite hybrid entanglement are considered. 

PACS numbers: 03.67.Mn 



I. INTRODUCTION 

In quantum-information science, the highly intriguing and 
nonclassical phenomenon of entanglement is at the heart of 
virtually all applications [1]. Being first considered by Ein- 
stein, Podolski, and Rosen, presenting the so-called EPR- 
paradox [2], entanglement has been in the focus of much 
research, especially in the recent past. However, there 
are stiU a lot of open questions remaining. When dealing 
with entanglement, one normally deals either with purely 
discrete-variable (DV) quantum states which live in a finite- 
dimensional Hilbert space, or with fully continuous-variable 
(CV), infinite-dimensional quantum systems. The toolbox for 
analysis is quite different for the two settings. While in the DV 
case, density matrices provide a complete and convenient rep- 
resentation, the investigation of CV states, due to their infinite 
dimensionality, is more subtle. However, at least in the Gaus- 
sian case, finitely many first and second moments are suffi- 
cient for a compact and complete representation [3-7]. 

Recently, so-called hybrid protocols have emerged which 
utilize both CV and DV resources at once [8-12]. These re- 
sources may include CV and DV states as well as CV and 
DV quantum operations and measurement techniques. It is 
worth noting that the term hybrid, which seems to be quite in 
vogue at the moment, is used in different contexts in quan- 
tum information. For example, there are proposals consid- 
ering hybrid quantum devices which combine elements from 
atomic and molecular physics as well as from quantum op- 
tics and also solid-state physics [13]. Furthermore, there is 
the notion of hybrid entanglement, referring to entanglement 
between different degrees of freedom, for example, the en- 
tanglement between spatial and polarization modes [14, 15]. 
However, in the present work, we use hybrid in the above 
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first-mentioned sense and define hybrid entanglement as the 
entanglement between a finite-dimensional, DV quantum sys- 
tem and an infinite-dimensional, CV quantum system. The 
prime example is an entangled state between an atomic spin 
and an electromagnetic mode. As already pointed out, the de- 
scription of CV and DV states would typically differ Hence, 
combining CV and DV quantum systems has its own charac- 
terizing and challenging subtleties. So, why would it be useful 
to consider such hybrid approaches? 

In CV quantum computation, Gaussian states as well as 
Gaussian transformations, such as beam splitting and squeez- 
ing, are used. However, to reach computational universal- 
ity just linear Gaussian elements are not sufficient [16]. At 
least one non-Gaussian component is necessary. Actually, any 
quantum computer utilizing only linear elements could be ef- 
ficiently simulated by a classical computer [17]. This single 
non-Gaussian element is the main challenge in CV quantum 
computation, as it is very difficult to efficiently realize such 
non-Gaussian transformations. In DV quantum computation 
with photons, the encoding of information takes place in a 
finite-dimensional subspace of the infinite-dimensional Fock 
space. Just as in the CV case, for deterministic processing, 
a nonlinear interaction is required to realize DV universal- 
ity [18]. But when truncating the Fock space, only single 
or few-photon states are left. The drawback of optical DV 
quantum computation then is that nonlinear interactions on 
the few-photon level are hard to achieve. Note that a well- 
known efficient protocol for universal DV computation with 
only linear optics is still probabilistic (or near-deterministic 
at the expense of complicated states entangled between suffi- 
ciently many photons) [19]. 

A way out of the problems of CV and DV quantum com- 
putation may be provided by those hybrid approaches. The 
scheme by Gottesman, Kitaev, and Preskill, which makes 
use of CV Gaussian states and transformations in combina- 
tion with DV photon number measurements, can be consid- 
ered one of the first hybrid protocols for quantum compu- 
tation [20]. So-called non-Gaussian phase states are created 
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from Gaussian two-mode squeezed states with the aid of pho- 
ton counting measurements. Additionally, the protocol can 
be considered as hybrid, as it employs the concept of en- 
coding logical DV qubits into CV harmonic oscillator modes 
(qumodes). 

However, there are other situations in which one may ben- 
efit from combining CV and DV techniques. For example, 
hybrid entangled states can be exploited for the generation of 
coherent-state superpositions within the framework of cavity 
QED [21]. Furthermore, there are quantum key distribution 
schemes which make use of hybrid entanglement [22-24]. Fi- 
nally, so-called qubus, i.e., quantum-bus-based, schemes have 
been developed for establishing entanglement between distant 
qubits. These also involve hybrid entanglement [25-27]. 

We can conclude that hybrid entanglement is a key ingredi- 
ent of various recent quantum-information protocols. In this 
paper, we perform a thorough classification of hybrid entan- 
gled states with the focus on bipartite hybrid entanglement. 
Nevertheless, also multipartite hybrid entangled states shall 
be briefly discussed at the end of the paper In addition to the 
derivation of a complete classification scheme, which distin- 
guishes between DV-like hybrid entanglement and true hybrid 
entanglement, also entanglement witnessing and quantifica- 
tion of hybrid entanglement will be discussed. 

The paper is organized as follows. In Sec. II, we derive 
the classification scheme for bipartite hybrid entangled states. 
First, the non-Gaussianity of hybrid entanglement is proved. 
Then we introduce an inverse Gram-Schmidt process, which 
is used for the derivation of the classification scheme. In Sec. 
Ill, examples for each class of hybrid entangled states are pre- 
sented. Section IV briefly discusses multipartite hybrid entan- 
glement, before we give a conclusion in Sec. V. Appendices 
A and B present auxiliary calculations. 

II. CLASSIFYING BIPARTITE HYBRID 
ENTANGLEMENT 

Let us start with a rigorous definition of bipartite hybrid 
entanglement. 

Definition 1. Any entangled bipartite state of the form 

N N 

n=l n=l 

d-l 

l¥«) = E C„m\m)^ \%,n,f , (1) 
m=0 

d-l 

c„„, eCVn,m; E kmnP^l, 

m=0 

for I < N < °° with generally nonorthogonal qumode state 
vectors |\|/„m)^, defined in the total Hilbert space 9{^^ — 
(E) 9-[£ with finite d >2, is called hybrid entangled. 

Since every bipartite hybrid entangled state can be written 
in such a pure-state decomposition, this definition is complete. 
Note that later, in Sec. II C, the rank parameter will deter- 
mine one of three possible classes of hybrid entangled states. 



A. Non-Gaussianity 

How can hybrid entangled states actually be described in a 
convenient way? Can their entanglement be quantified? For 
an overall infinite-dimensional Hilbert space (g) !J-(^, us- 
ing the standard finite-dimensional techniques, density matri- 
ces can no longer be employed. So, the use of CV methods 
can be attempted. However, in CV entanglement theory, the 
only conveniently representable states are the Gaussian ones. 
In the non-Gaussian CV regime, for instance, exact entangle- 
ment quantification is, in general, hard to achieve. It is there- 
fore useful to ask whether hybrid entangled Gaussian states 
exist. 

Lemma 1. Any single-partite d -dimensional quantum state 
with finite d and d >2 is non-Gaussian. 

For the proof of Lemma 1, see Appendix A. 

Tlieorem 1. Any bipartite hybrid entangled or classically 
correlated state is non-Gaussian [46]. 

Proof. If a multipartite quantum state is Gaussian, all its 
subsystems will be Gaussian. So, we consider a state of the 
form (1) and trace out the CV subsystem. What is left is a d- 
dimensional single-partite system with finite d, which can be 
described in its Fock basis. It is denoted by p. Due to Lemma 
1, p is Gaussian if and only if p = |0) (0|, which is a pure 
state. However, for any entangled or even classically corre- 
lated state, the reduced state cannot be pure [28]. Hence, for 
Gaussian p, the overall system cannot be entangled or classi- 
cally correlated. Therefore, every bipartite hybrid entangled 
or classically correlated state is non-Gaussian. ■ 

We have found that there are no hybrid entangled Gaussian 
states. Hence, CV Gaussian tools are inappropriate for the 
description of hybrid entangled states. 

The proof of Theorem 1 basically relies on Lemma 1, which 
states that any finite-dimensional DV state with dimension > 2 
is non-Gaussian. It is straightforwardly generalized to multi- 
partite systems. Arguing that for Gaussianity all subsystems 
have to be Gaussian, it is sufficient for the non-Gaussianity 
of any multipartite system which also possesses DV con- 
stituent(s) that at least one DV subsystem is of dimension 
> 2. However, if there shall be entanglement between the 
DV subsystem and the rest, the state will necessarily be non- 
Gaussian, since for entanglement dimension, > 2 is required. 
So, any multipartite quantum state which involves entangled 
DV subsystems is non-Gaussian. Even more generally, only 
1- or infinite-dimensional systems can be Gaussian. This is, of 
course, not surprising; however, it is worth pointing out that 
Theorem 1 is not trivial. Since the overall Hilbert space of 
hybrid entangled systems is indeed infinite-dimensional, it is 
not a priori clear that Gaussian hybrid entangled states do not 
exist. 

We have shown that hybrid entangled states belong to 
the non-Gaussian, infinite-dimensional Hilbert-space regime, 
which is not easy to deal with, as we already know from con- 
ventional CV entanglement theory. The states can be neither 
described by proper density matrices nor by covariance matri- 
ces. Phase-space representations are also not so convenient. 
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since one of the subsystems is DV. The only known quasi- 
probability-distribution which may in some cases make di- 
rect statements about the separability properties of the state 
is the Glauber-Sudarshan P representation [29, 30]. How- 
ever, it can be easily shown that this function is totally ir- 
regular for these highly nonclassical hybrid entangled states. 
It may still be possible to construct entanglement witnesses, 
but entanglement quantification appears, in general, hard. Ac- 
tually, there is a way out of this dilemma: For some hybrid 
entangled states, the unique Hilbert-space structure can be 
exploited in such a way that the states can nevertheless be 
described by density matrices. These states are effectively 
finite-dimensional and are therefore called DV-like hybrid en- 
tangled. This gives rise to a classification scheme of hybrid 
entangled states. However, for a rigorous analysis we first 
need to introduce a slightly modified version of the so-called 
Gram-Schmidt process. 



B. Inverse Gram-Schmidt process 

Consider again a general hybrid entangled state of the form 
(1). Depending on the number of mix terms and the dimen- 
sion of the DV subsystem d, there can be maximally N x d 
linearly independent CV qumode states |\|/,„„)^ in p"^^. The 
dimension d is always finite due to the definition of hybrid 
entanglement. However, may be either finite or infinite and, 
hence, the number of linearly independent qumode states is 
either finite or infinite. Furthermore, if = 1, the state is 
pure. 



If the number N x d of linearly independent CV qumode 
states is finite, they only span an (A^ x t/) -dimensional sub- 
space y^Nxd of the initially infinite-dimensional Hilbert space 
y-L. Then, the Gram-Schmidt process can be employed to ex- 
press the qumode states in an orthonormal basis of this finite- 
dimensional subspace. In this case, the state becomes effec- 
tively DV and all the methods from DV entanglement theory 
can be applied. 

The Gram-Schmidt procedure is a method for orthonormal- 
izing a finite, linearly independent set of vectors in an inner 
product space [31]. For a linearly independent set of vectors 
{|\|/,) : i— !,...,«} (since the process is to be exploited in the 
framework of Hilbert spaces, the inner product space is a pri- 
ori assumed to be a Hilbert space), a set of pairwise orthonor- 
mal vectors {|e,) : i — !,...,«} spanning the same subspace 
as { |\(/,) : / = 1 , . . . , n} is given by 



k'l) = l¥i) , 



14) = l¥2)-(ei|¥2)ki 



7(4 



(2) 



W„) = - L {e,Wn) let) , \e„) = -Jfi^ . 

Making use of this, any finite set of linearly independent 



qumode states can be expressed in an orthonormal basis { |e,) : 
i — I,... ,n}. However, more specifically, Eqs. (2) only deter- 
mine how to express the new orthonormal basis in terms of 
the old nonorthonormal one. What is actually required is the 
inverse expression. To express the qumode states in terms of 
{|e,) : i — 1, ... ,n} we try the following approach, which can 
be considered an inverse Gram-Schmidt process, i.e., a modi- 
fied version of the original Gram-Schmidt process. 

Theorem 2. n normalized, generally nonorthogonal, lin- 
early independent states {|\(/,) : / = 1, . . . ,n; < | | < 
1 V;, j} can always be expressed as |\(/,) — L;=i <^ii kj). where 
{ : ; = 1 , . . . , n} forms an orthonormal basis of the space 
spanned by {|v|/,)}, andatj G C. 

Proof. Writing out |\(/,) — Y!i=i'^ij\^i)^ Theorem 2 states 
that the {IV;)} can be always written as 

l¥i> =flii ki) , 

|\|/2) =021 kl)+«22k2) , 

|\(/3) = fl31 kl) +"32 k2) +fl33 k3) , 



(3) 



|x|/„) =fl„i \ei) 



1=1 



For the proof, it has to be shown that (1) |\(/,) = L}=i '^ij k/) 
corresponds to a valid basis transformation and (2) it actually 
performs the right mapping. 
(1) Write the transformation as 



Wi)=Y.Tij\ej) 



(4) 



with the transformation matrix 
fan 

021 fl22 



(5) 



1 ^/i2 ' ' ' a^jYij 



Due to the normalization of the initial and the new vectors 



1 1"'; I 



1 is known, and due to the linear independence 



of the also an ^ OV/ is necessary. 

det[r]-nLifln^O. 
r is invertible. 
=^ r is a valid basis transformation. 

(2) To show that the lower triangular structure of the basis 
transformation T in combination with the orthonormal basis 
{|e,)} is sufficient to actually express the {|v(/()} accurately in 
terms of {k/)}' it is demonstrated that the ^li^^tii parameters 
fl,y can be chosen such that all overlaps (\(/, |\(/,) are preserved 
when the transformation is applied. 

On the one hand, there are n^ such overlaps in total and 
"^"2^^^ ones with potentially differing absolute values. On the 
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other hand, there are complex parameters fl,y. From 

the structure of the basis transformation and the fact that ajj 
are complex, it is clear that if aij can be chosen such that the 

^^^^y^ -element set of : / < j} can be preserved, also 

the rest of the overlaps are preserved, since they are only com- 
plex conjugates of the former. Hence, it already becomes rea- 
sonable that the a,y can be chosen appropriately. 

However, a proper proof is performed by induction in n: 
Inductive Basis: n = 1. There is only one overlap to be 
preserved: 



(6) 



Hence, choose an — I, which preserves the overlap 

Inductive Step: Assume aij have been calculated for i <n — 
1 such that all overlaps : i,j <n~ 1} are preserved. 

We show that then also all a„j can be chosen such that the 
overlaps {(\(/,|\(/„) : / = 1, . . . ,n} are preserved. The complex 
conjugated overlaps follow automatically as argued before. 

Applying the basis transformation, the overlaps {(\(/, |x|/„) : 
; = l,...,n— 1} are 



(¥l|V«) =flllflnl , 
(¥2|¥«> = fl21«nl +a22«n2 , 

(¥«-l l¥«) ~ ^*i-l,l^nl + • • ■+^*i-l,n-l^n.,n-l ■ 



(7) 



As atj have been calculated for / < n — 1 due to the inductive 
hypothesis, this is just a system of linear equations, which can 
be written as an augmented matrix (using a\^ = 1): 



/ (¥il¥"> 
(¥2|¥«) 

V(¥«-il¥«) 



1 


.. 











"21 


«22 





(8) 



■'n-lA "n-1,2 



Since a,, ^ OV/, this system of equations is exactly solv- 
able. Hence, {a„j : j — I, . . . ,n — 1} can be chosen such that 
{(\|/,|v|/„) : ! = 1, . . . ,n— 1} are preserved. Therefore, there are 
only one free parameter a„„ and one overlap (x|/„|x|/„) to be 
preserved left: 



(¥«l¥«) = 1 = E 

7=1 



n-l 



njl 



(9) 



From T!j=i Wnj?' = 1 ^nd a„„ ^ 0, which is already known, 
Lj=i Wnjl^ < 1 follows, and hence «„„ can be chosen as 



\ 



7=1 



(10) 



In the end, also (x|/„|x|/„) can be preserved. Therefore, the the- 
orem is valid for n under the assumption of validity for n — 1. 
As a conclusion, with the inductive basis it is valid for all n. 



The proof has been presented in such a great detail, because 
it is constructive and hence also sets out how to actually com- 
pute {aij} for a given set of qumode states. Recalling Eqs. 
(3) and (8), ajj can be calculated successively one after an- 
other by considering successive overlaps. A parameter an is 
directly obtained from the overlap |\(/,). Then, a,2 follows 
from (i|/2|\|/,) together with the known an. Likewise, 0,3 is 
calculated from (\(/3 an and oq- For the other parameters 
just go on like this. Hence, the inverse Gram-Schmidt process 
can be efficiently implemented and computed. 

As an example consider the normalized qumode states 
:/= 1,2.3} with overlaps 



(¥il¥2) =ci, 
(¥il¥3) =C2, 

(¥2|¥3) =C3. 



(11) 



They can be expressed in an orthonormal basis {|e,) : / 
1,2,3, {ei\ej) = dij} as 



l¥ 



1 = ei 



l¥2) =ci|ei) + ^l-|ci|2|e2) 
l¥3) =C2|el) + ^==^=|e2) 



(12) 



|2 



'i-k2 



|C3-C,C2 
l-hP 



k3). 



C. Classification scheme 

So, once again, consider a state of the form (1). For N ^ I 
the state is pure. Then it contains only d linearly indepen- 
dent qumode states, which span a tZ-dimensional subspace 
Therefore, with aid of the inverse Gram-Schmidt pro- 
cess these qumode states can be expressed in an orthonormal 
basis. Then, density matrices can be employed for the descrip- 
tion of the overall state, and also a Schmidt decomposition can 
be performed or pure-state measures such as the entropy of 
entanglement can be calculated [28, 32]. The state is then ef- 
fectively DV. For 1 <N <°°, the state is mixed. Nevertheless, 
it possesses a finite number of Nxd qumode states, which can 
be again cast in an orthonormal basis. Therefore, also states 
of this kind are effectively DV and the density matrix formal- 
ism can be exploited. However, such states are no longer pure 
and neither pure-state measures nor a Schmidt decomposition 
can be applied. Finally, there is the case N Here, N = °° 
refers to those states which can be expressed with infinite 
only. These states hold an infinite number of qumode states, 
which has the effect that the Gram-Schmidt process cannot be 
applied anymore. Hence, they are not describable by density 
matrices and therefore no longer effectively DV. Nevertheless, 
one subsystem does remain DV. These are the states which we 
call truly hybrid entangled. 

Summing up, for bipartite hybrid entangled states in a pure- 
state decomposition with denoting the number of mix terms 
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in the convex combination of pure-state projectors, there is the 
following classification scheme: 

9 N — 1 : Pure hybrid entangled states. 

- Supported by a finite-dimensional subspace. 

- DV-like entanglement. 

- Schmidt decomposition applicable. 

- DV pure-state measures applicable. 

• 1 < TV < cxD : Mixed hybrid entangled states. 

- Supported by a finite-dimensional subspace. 

- DV-like entanglement. 

- DV mixed-state measures applicable. 

• N — oo : Mixed hybrid entangled states. 

- No support by a finite-dimensional subspace. 

- => True hybrid entanglement. 

- In general no exact measures directly applicable. 

- CV entanglement witnesses adaptable. 

It should be pointed out that the possibility of applying DV 
methods on such a wide class of hybrid entangled states 
is quite remarkable. An initially non-Gaussian infinite- 
dimensional quantum state, which seems rather awkward at 
first sight, can finally be conveniently described in terms of 
density matrices, and in the pure-state case even a Schmidt de- 
composition can be performed. Nevertheless, there is also the 
class of states which stay truly hybrid entangled and cannot be 
transformed using the Gram-Schmidt process. As mentioned 
earlier, in this infinite-dimensional, non-Gaussian regime ex- 
act entanglement quantification appears to be hard. However, 
entanglement witnesses for detecting true hybrid entangle- 
ment can be adapted from CV entanglement theory. 



III. EXAMPLES AND APPLICATIONS 

Now, we shall present an example for each class of hybrid 
entanglement. We shall also point out when these examples 
correspond to manifestations of entangled states encountered 
in quantum-information protocols and apphcations. 



(1) An inverse Gram-Schmidt process with respect to sub- 
system B yields 



with X 

\AB 



\a)'' =x\eo)'' + ^/Y^\e,)\ 
|-a>^=x|eo>^-xVl-x2|ei)^ 



1 I |2 

e~i\^\ . This corresponds to Eq. (12). Hence 



(14) 
(15) 

(16) 



IV)^^ =;7f (l^o)^ l^o)'' +x \e,t |eo)« + VI^^\e,t \e,Y 



-x\e2)^ \eQ)^ -x^\/\^ 



-x2-x4+x6|e2)^|e2)^ 



(17) 

which is the effective DV form of the state. 

(2) The pure-state entropy of entanglement of the state can 
be calculated, which is shown in Fig. 1. For (X^ °° the state 
becomes maximally entangled (a real). 

Entropy of Entanglement 

Towards Maximal Entanglement 




Amplitude a 



FIG. 1: (Color online) Entropy of entanglement of the state of 
Eq. (13). For a — > oo the state becomes maximally entangled. Note 
that the entropy of entanglement has been calculated here in qubit 
entanglement units (ebits: log = logi)- 



(3) Finally, set x = j, which corresponds to a = \/21n2 ; 
1.18, and calculate the Schmidt decomposition: 



/\A| i.B 



Wf'' =Q.16\eX\e',r +0.56\e\)^\e' 



-0.33|4)^|4)^ 



(18) 



We can also think of other examples for pure bipartite hybrid 
entangled states. States of the form 



A. Pure qutrit-qumode entanglement 

As an example of pure bipartite hybrid entanglement an en- 
tangled state of a qutrit and a qumode system is considered: 

l¥)^^ - ^ (ko)^ |vac)^ + ki)'^ |a)^ + \e2f |-a)^) . (13) 



|v)^^ = i=(ko)^|Vo)" + h)^l¥i>") 



(19) 



are pure bipartite hybrid entangled and relevant for cat- 
state engineering [8, 21], hybrid quantum communication via 
qubus approaches [25, 27], and some quantum key distribu- 
tion schemes [22, 23] [later, in Sec. IIIC, we obtain such a 
state in Eq. (45) as a special case of Eq. (26) with (n,/,) ~ 
and r] = 1]. 



6 



B. Mixed qubit-qumode entanglement 

Now, we present an example for the class of mixed, but 
effectively DV hybrid entangled states. Consider the state 



\AB 



\AB , 



>^|vflc)^ + |ei)^|±a)^ 



(20) 



which contains three qumode states \vac)^ and |±a)^ with 
a e R (see Fig. 2 for a visualization). 



I0> 
^vac> — 



ll> 

— |-0f>- 



ll> 



1-p 



FIG. 2: (Color online) Visualization of the hybrid entangled state 
p'** of Eq. (20) in phase space. The blue (right) region corresponds 
to the first pure state in the convex combination with probability p, 
while the green (left) region represents the pure state obtained with 
probability 1 — p. The additional |0) and 1 1) vectors denote the qubit 
states associated with the qumode states \vac) and |±a). 



We perform a Gram-Schmidt process and obtain a qubit- 
qutrit entangled state in 9-C2 ® ^f- The pure states in the 
convex combination of (20), after the inverse Gram-Schmidt 
process in an orthonormal basis, look like 



1^- 



\AB 



\AB 



1 / 



V2V 

u 



( ko)"^ |eo>^ +x\ei)'^ |eo>^ + V 1 -x^ ki)'^ kl>^ 

(21) 



ko)"^ \eo)^ +x\ei)^ \eo)^ -x^\/ 1 -x^ {eiY 



'1} 



(22) 



with X = e^zl"! . Entanglement quantification can be per- 
formed with the logarithmic negativity En (see Fig. 3) 
[33, 34]. For general bipartite mixed states with higher di- 
mension than 2 X 2, it is basically the only known calculable 
entanglement monotone, and it is also the most important one 
for general mixed, effectively DV hybrid entangled states. It 
can be seen that for any given a the smallest entanglement is 
obtained for maximal mixing p = \. Furthermore, the higher 
a, the more entangled the state is. Finally, for = V 1 and 
a 00 the state approaches maximal entanglement. 




0.5 

Probability p 



Amplitude a 



FIG. 3: (Color online) Logarithmic negativity of the effective qubit- 
qutrit hybrid entangled state p"^^, Eq. (20), as a function of proba- 
bility p and amplitude a. For any given a maximal mixing p = 2 
yields the smallest entanglement. Furthermore, higher a also results 
in greater entanglement. For p = V 1 and ^ the state becomes 
maximally entangled (a e R). 



We may also consider other mixed bipartite hybrid entan- 
gled states which are even effectively DV-like qubit-qubit en- 
tangled. Then the concurrence [35] is a better choice for en- 
tanglement quantification. In Sec. Ill C, such a state is given 
in Eq. (43) as a special case of Eq. (26) for («(/,) = 0. Its 
concurrence is plotted in Fig. 8. 



C. True hybrid entanglement 

Finally, the third class of bipartite hybrid entangled states 
is considered, the truly hybrid entangled states. Recall their 
definition corresponding to Eq. (1) with N = °°: 



P'^^ = L P" I^")ab (v« I , /?„ > V « , 



m=0 



n=l 
ln=0 



(23) 

As can be seen from these equations, truly hybrid entangled 
states possess an infinite number of qumode states |x|/„m)^- 
Therefore, the Gram-Schmidt process ceases to work and 
the states stay in a Hilbert space of the form ^H:f (g) H^. 
Hence, they are not effectively DV, but instead really com- 
bined DV (E) CV states and therefore truly hybrid. Unfortu- 
nately, this true "hybridness" has the effect that the states live 
in an overall infinite-dimensional Hilbert space in the non- 
Gaussian regime. This makes exact entanglement quantifica- 
tion difficult. Now, we present an example which is relevant 
for some quantum information protocols. 
Consider the state 



AB 



\af + \l)^\-af), (24) 
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which plays a crucial role in cat-state engineering [21], hy- 
brid quantum communication via qubus approaches [25, 27], 
and some quantum key distribution schemes [22, 23]. It is 
transmitted through a one-sided thermal photon noise chan- 
nel, where the noise only affects subsystem B. Writing this 
out, 

p'^« = (l'^®$L„,„,)|vr(v| (25) 

is the state to be investigated. The channel is modeled by a 
beam splitter coupling subsystem B to the environment, which 
is in a thermal state. Afterward the environment is traced out 
(see Fig. 4). 




FIG. 4: (Color online) Modeling of the photon noise channel. The 
environment mode is in a thermal state and coupled via a beam split- 
ter to the input state to be transmitted. Subsequently the environment 
mode is traced out and the decohered output state is obtained. 

For the output, 

p'^'=(i-'®$L™„)ivr(¥i 

^ n=0 A:,/=0 




is obtained, where a and a ' are the mode operators of subsys- 



tem B; p„ denotes the thermal photon distribution ^^_^_j."'* j^„+i 
of the environmental thermal state with mean thermal photon 
number («,/,), and/„/.(r|) andA^|^(r|) are defined as 

UM) ^ Q ^/n"-'(-y^^)^ (27) 

Klii^) ■■= (yr^Pi«"""«^""'iv^«) , (28) 

with the beam splitter transmissivity r| (for a detailed cal- 
culation of this result, see Appendix B). It can be inferred 
from Eq. (26) that the state p'"^^ clearly is truly hybrid en- 
tangled, as it contains an infinite number of qumode states 

a- A: 

{a' |±-yf|a)g : A: = 0, 1, . . • ,°°}. Hence, true hybrid qubit- 
qumode entanglement is obtained. Furthermore, the form of 
the output state illustrates the effect of the thermal photon 
noise channel in a concrete way: On the one hand, the damp- 
ing effect due to the beam splitter is clearly visible in terms 
of yfi in the states d' \±^Jf[a). On the other hand, there 
is not only damping, but also thermal photon noise, which 
becomes manifest in the creation operators in the states 
fl' |±-yr|a). Thermal photons "leak into the system" and are 
"created" in the damped coherent states. Finally, it is clearly 
visible how each term of p^,g„„„; = £^=0 P'n 1")^ ("I results in 
the creation of at most n noise photons in the coherent states. 
However, as descriptive as Eq. (26) is, it is inapplicable for 
further calculations. For example, for entanglement witness- 
ing, moments may have to be computed [37, 38]. Unfortu- 
nately, in this case, making use of the state as written in Eq. 
(26), intractable infinite sums are obtained whose convergence 
behavior is impossible to be worked out exactly. Of course, 
truncation at some value n could be performed, which would 
certainly result in very accurate outcomes provided that value 
is large enough [36]. However, such a procedure is opposed 
to the actual intention of analyzing true hybrid entanglement, 
since a truncated state is not truly hybrid entangled anymore. 
This is a point which makes the investigation of true hybrid 
entanglement particularly challenging. Infinite sums or inte- 
grals emerge, which have to be calculated exactly. 

However, the transmitted state can be also written as 



+ |i)^(i|®l-ViTa-v/^^y)^(-\/iT«-v/^^Yl 

+ A-««^(ii) 1 1)^ (oi ® \-^a ^r^yf (v^oc - yr^Yi) , 
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with 

A«pY(n) 



1 , 



exp 



(30) 



- ( /r^p* + Vnf ) ( /T^a + Vny) 



(see Appendix B for the calculation). This form of the state 
p/AB jjQj. insightful as that of Eq. (26). However, it 
is mathematically much more convenient. Instead of infinite 
sums, in this form integrals occur which make the calculation 
of moments straightforward. 

Now we want to derive entanglement witnesses for true hy- 
brid entanglement. For this purpose, we exploit the determi- 
nant (now a and are the mode operators corresponding to 
system A and b and b'' those of system B) 



1 (a 

{a) 



j^b) 



{a^ a) 
a' ah' 



(a'ab 



ib^b) 



which Shchukin and Vogel (S V) proposed and utilized in their 
work [37, 38]. However, we have to slightly adapt the SV ap- 
proach to our case, since, instead of a fully CV system, we 
deal with a hybrid quantum system of a DV and a CV subsys- 
tem. There are two ways for this adaption. 

The first way is to simply interpret the DV subsystem as 
living in a subspace of an infinite-dimensional Hilbert space. 
The DV qudit is interpreted as a CV system supported by !J-L 
and being encoded in a Fock basis. However, it only makes 
use of a finite number of the basis vectors. Then the SV crite- 
ria can be applied just as usual. 

The second approach is not to adapt the state, but to adapt 
the criteria. Assume the system to which the operators d and 
fl^ belong is c/-dimensional. Then the orthonormal Hilbert- 
space basis vectors \m) can be written as column vectors with 
a "1" in row m with m = 0, . . . ,(i — 1: 



= (Oo, ■ ■ ■ ,Om-l, 1,„,0„,+ 1, . . . ,Orf_l)' 



(32) 



With this notation the new qudit mode operators dii and aJi can 
be defined as proper d x d matrices: 



/O Vl 




Vo 



\ 



/ 
... 0\ 



^ 
: V2 



VcT^ 0/ 



They have the properties 



ad \n) 



« In — 1) 



|n) = (l-§«,^-i)V^H^|n+l): 



i^dY = 0, 

(4)' = o. 

Furthermore, the commutator becomes 



(35) 
(36) 
(37) 
(38) 



[dd,aj = 



/I 



td- 



-{d-l)J 



(39) 



dxd 



dxd 



Since any t/-dimensional qudit \\\t)^i can be written as \\\t)^i ~ 
^^(arf)|0)j for an appropriate operator function g — g{aii) 
(31) and also the projection operator on |0)^/(0| can still be ex- 
pressed as : e^"d"'' :, where : . . . : denotes normal ordering 
[37], the SV criteria can be also derived for hybrid systems 
with aid of these new operators. Structurally the same criteria 
are obtained, except that for certain 11,12, jiji the moments 
M.jip) 



[a J a^cij aj t>' b b^ b )- are zero. Due to the 
properties (37) and (38) any combination of /i,;2,ii,i2 such 
that {cidY or {a^fY with k>d occurs in the moments, nulli- 
fies them. In the end, in the matrix of moments the rows and 
columns corresponding to these combinations of 11,12, jiji 
are simply missing. 

When applying the SV criteria it hardly makes a difference 
whether the first or the second approach to the adaption of 
the criteria is chosen. Only for moments involving terms like 
(fl^)*^(flj;)' with A:, Z e INJ the two approaches may yield differ- 
ent results. However, for the determinants used in this paper, 
both lead to the same result. 

For the calculation of the SV determinant s in Eq. (31) for 
the state in Eq. (29), 



ipMB(a,ri, 



l-Tl 



{nth) 



-4 ah 



(40) 

is obtained. A graphical illustration for r| = | is shown in Fig. 
5. It can be observed that there clearly is a parameter region 
in which entanglement can be detected. Also note the trade- 
(33) behavior and the optimal a"'" which corresponds to the 

most robust state p'^o,,, regarding entanglement witnessing for 
fixed T| = I and varying mean thermal photon number («(/,). 
The origin of this trade-off is based on the trade-off between 
too little initial entanglement before the channel for low a and 
very fragile entanglement for too high a (a G R). 

Furthermore, it is worth pointing out that the witnessed en- 

(34) 

tanglement in this case is actually true hybrid entanglement. It 
can be concluded that the SV determinants provide a suitable 
tool for the detection of true hybrid entanglement. Note that 
we exploited determinants involving moments of 4''' order 
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FIG. 5; (Color online) The upper graph displays the SV determinant 
i'p/As(a,Ti = 5, (n,ft)) for the state p''*^ with ^ = ^- Without loss of 
generality a e R has been assumed. On the lower diagram the two 
regimes are plotted. There can be clearly identified a region, where 
the determinant is below zero. Hence, true hybrid entanglement can 
be witnessed. Furthermore, there is a trade-off behavior which results 
in the existence of an optimal a"'". This a"'" corresponds to the most 
robust state p'^?,„ regarding entanglement witnessing for fixed 11=1 
and varying mean thermal photon number («(/,). 



There are also other approaches for witnessing entanglement 
in similar states, which involve lower-order moments [23, 39]. 
However, often it is necessary to consider moments of higher 
order to detect entanglement. For example, the detection of 
entanglement in the pure state iA(^(a)(|a, a) — |— a, — a)) re- 
quires moments of 4''' order [37]. 

When setting 5a/as (oc, r| , {n,h) ) < and solving this inequal- 



ity for {n,h) 



(nth) < 



4ri|a|' 



(1-Tl)(2e^ 



4|a|' 



1) 



(41) 



is obtained [47]. For parameters (a,r|, (n,/,)) satisfying this 
inequality, entanglement is detected. Furthermore, the in- 
equality can be used to define a surface. The parameters 
(a, Tj, («(/,))(.„? for which entanglement is verified lie below 
this surface (see Fig. 6). However, it is rather cumbersome to 
read off exact parameters in such a 3D plot. Hence, regions of 
successful entanglement detection are plotted in Fig. 7 for dif- 
ferent values of r|. As expected, the higher the transmissivity 
r| the greater the parameter regions of entanglement detection. 
Furthermore, we can infer from Eq. (41) that the optimal a"'" 



Entanglement Verified For Values Below Surface 




Thermal Photons <n> 



Transmission jy 



1.0 

1.5 Amplitude a 



FIG. 6: (Color online) Surface defined by Eq. (41) (a £ R). For 
parameter triples (a,r|, (n,/,)) lying below it, entanglement is wit- 
nessed. Once again, the trade-off behavior can be recognized and an 
optimal a"^" for which entanglement can be detected in the presence 
of the strongest possible noise. 
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FIG. 7: (Color online) Regions of entanglement detection for differ- 
ent values of transmissivity r|. Note that the optimal a"'" has a fixed 
value (a e R). 



does not depend on rj. This is also recognizable in Fig. 7. We 
find a"'" « 0.44. It is actually quite remarkable that the op- 
timal amplitude a"'^' regarding entanglement witnessing with 
the SV determinant s does not depend on the channel param- 
eters at all; a"'" is determined solely by the choice of the SV 
determinant. 

Compare this to Fig. 8, which shows the concurrence [35] 
of the state when setting (n,/,) = 0, which is given by 



C(p'^^) =^ - e-4ill«l- ( + 3e-4(i-^)l«|- 



-4(l-r|)|a|- 



(42) 



There, the optimal a"'''(r|) does depend on the transmissivity 
r|: The greater the transmissivity, the higher a"'"(r|). This 
is quite remarkable. It can be inferred that the ability to de- 
tect entanglement in p''*^, depending only on the choice of 
the SV determinant, does not behave in the same way as the 
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Concurrence 
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FIG. 8: (Color online) Concurrence of the state p''** for (n,/,) = 0, 
as a function of transmissivity r| and amplitude a. On the one hand, 
the higher a the greater the initial entanglement, but also the more 
sensitive the state is to photon loss. On the other hand, for low a 
there is only little initial entanglement. However, the state is more 
robust against losses. Hence, there is a trade-off behavior, and an op- 
timal a!^" depending on T| exists for which the output entanglement 
is maximal. 



entanglement of p'^^ itself, which of course depends on the 
thermal channel's parameters r| and («(/,). 

Note that in this case, the state of Eqs. (24) and (25) actu- 
ally becomes an effectively DV hybrid entangled state, when 
setting the mean thermal photon number of the channel to 
zero. The state is then subject only to amplitude damping, 
resulting in 



1 



+ \l)^{l\®\-^af{-V^a\ 



,-2(l-ri)|a|' 



|0)^(l|®|^a)^-V^a| 



(43) 



-e-2(i-^)l«r ID 



which is mixed but contains only two qumode states 
|±y^a)^. Hence, we can perform an inverse Gram-Schmidt 
process and describe the state by a proper effective qubit-qubit 
density matrix: 



^ {n,i,)=0 ■ 







1 Xk 


he A.2 XVl^ 

\kVT^ X^/l^ 1-A.2 J 



(44) 



withK:=e-2(i-ri)|aF andA.:= (-0Ta| Vrfa) = e-^nlap^ The 
concurrence of this state is given in Eq. (42) and plotted in 
Fig. 8 for a e R. 

If, besides setting («,/,) = 0, we also set the transmissivity 
of the channel to one, the channel is completely canceled and 
we obtain a pure bipartite hybrid entangled state: 



l¥) 



AB 



U|0>^|a)^ + |l)^ha)^). (45) 



V2 



Hence, in this physical example, we can access all three 
classes of bipartite hybrid entanglement and even switch be- 
tween them using the transmissivity r] and the mean thermal 
photon number (nr/i)- 



D. Yet another truly hybrid entangled state 

In order to demonstrate that we can easily construct other 
quantum states featuring true hybrid entanglement, we shall 
present another truly hybrid entangled state and its entangle- 
ment detection: 



.AB 



Pn 



)^{\Qt\Vna) 



V2 



<x < 1 . 



/na) 



(46) 



As this state contains an infinite set of qumode states 
{li-y/noc)"^ : M = 1,2, . . . ,oo}, it is clearly truly hybrid entan- 
gled, exhibiting (true) hybrid entanglement between a qumode 
and a qubit. 

Now, we shall witness entanglement in this state with aid 
of the SV determinant s of Eq. (31). For the calculation of the 
moments, the identities 



1 



n=l 



E«J 



X 



are exploited, and we obtain 



(47) 
(48) 



s{x, oc) = - 



2a^ /a(l— x) 

1 —X \ X 



EV^(xe-2«^)" 



e-2«'(i-x)N^ ^a(l-.=c 
1 — xe^2a2 

a(l -x) ^ 



E^^(xe-2«^)« 



=1 



n=l 



a /e 
1-xV l_jt 



1-x) 



It is clear that 



for < X < 1. Hence, 



n=l 



xe 



< 



1 — xe 



xe 



< 



xe 



-2a- 



(l-xe-2«')2' 



(49) 
(50) 

(51) 
(52) 
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Inserting the lower bounds into the sums of Eq. (49) yields 



IV. MULTIPARTITE HYBRID ENTANGLEMENT 



s'{x,a) 



2x 



1-x 



-4a- 



1 



-1 -xe-2aV " V ' 1 -xy 

(53) 

Since s'{x,a) > s(jt:,a)Va e R,x e]0, 1[, from s'{x,a) < 
follows s{x,a) < 0. Therefore, s\x,a) < is a sufficient 
criterion for entanglement detection. It is plotted in Fig. 9. 
The graphics show that entanglement can be verified for suffi- 




1.0 

1.5 Amplitude a 



Probability x 
0.8 



In this section, investigations regarding multipartite hybrid 
entanglement are presented. In general, -partite hybrid en- 
tangled quantum systems live in Hilbert spaces of the form 
y-Ci^ (g) . . . (g) y-^i^, where some di are finite and some infi- 
nite. Considering tripartite hybrid entanglement there are two 
cases: Either the Hilbert space looks like iJ-Q^ (g) y-Q^ (g !J-L or 
y-Qi <E) y-L (E) Hx,, with finite di,d2- 

Consider the first case, where only one subsystem is CV. A 



general hybrid entangled pure state in 
be defined as 



) !}{^ may 



d 1 ,ch 



where |vc,j)*' represent some qumode states. As there are at 
most di X d2 such qumode states, a Gram-Schmidt process 
can be executed to write the state as a multipartite, effectively 
DV hybrid entangled state. Therefore, if all are linearly 

independent the state effectively lives in a Hilbert space of the 

In the second case, where the initial Hilbert space looks like 
!}{^ (E) a general hybrid entangled pure state is 



l¥) 



ABC 



d.oo 



\eir\<^uf\^L 



(55) 
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FIG. 9: (Color online) The upper graph shows s'(x,a) which is de- 
rived from the original SV determinant s{x,a). The yellow plane 
denotes zero. The lower diagram presents the cutting line between 
s'{x, a) and the zero plane. It displays the witnessing region. We also 
want to point out that for a = the state is of course not entangled 
and s'{x,a = 0) becomes zero. 



ciently small x and a. This can be understood in the following 
way. On the one hand, small x corresponds to little mixing be- 
tween the pure states jvi/,,)"^^, which themselves are highly en- 
tangled, depending on the amplitude y/na. For pure states of 
the form \\\i„) , entanglement witnessing can be performed 
perfectly with the SV determinant s. On the other hand, we 
have aheady seen in previous calculations (see Figs. 5 and 7) 
that entanglement detection via s fails for large a. 

To conclude, a second example for a truly hybrid entan- 
gled state has been presented whose entanglement can again 
be verified with aid of the SV criteria. 



where both and |v(/,;)'' represent qumode states, and the 
Cjj are chosen such that (v(/|v(/) = 1. Obviously, in this case 
an infinite number of qumode states is present for an infinite 
number of Cij ^ 0. No Gram-Schmidt process can be per- 
formed. Hence, these states show multipartite true hybrid en- 
tanglement. 

To conclude, in the multipartite regime, already for pure 
states there are characteristic differences between the differ- 
ent possible configurations. Furthermore, in contrast to the 
bipartite setting, there is also pure true hybrid entanglement. 
When dealing with two parties, true hybrid entanglement only 
occurs for specific types of mixed states, as described in the 
previous sections. This can now be generalized. True hybrid 
entanglement can be obtained in two ways: Either the sys- 
tem is mixed with an infinite number of mix terms; then it 
is sufficient that only one subsystem is CV. Or two or more 
subsystems are CV; then the state is not even required to be 
mixed. However, every hybrid entangled mixed state with a 
finite number of mix terms, which contains only one CV sub- 
system is effectively DV. An A^-partite mixed state of this type 
in (g) . . . (g) •!^^_| ® Hx, can be always described in a Hilbert 



space of the form O-Q^ 



®Hl^_^®^z, where 

N-l 

MY{di, 

i=l 



(56) 



and M is the number of mix terms. 

As an example, we want to discuss an explicit tripartite hy- 
brid entangled state, supported by a Hilbert space of the form 



12 



Hi<^^<^ with finite d = 2: 
l¥)^''^ - ^ (ko>^ Ic^o)'' l¥o)^ + ki)^ |c^i)^ l¥i>^) . (57) 

In this example, the special case occurs that there does not 
exist an infinite number of 7^ 0. In terms of Eq. (55), c,y — 
■Jj for / = 1,2 and 7 = 1, and c,,- = V j > I. Hence, there is 
no true hybrid entanglement and the state can be described by 
DV methods. 

Defining Qj, := ((|)o|(|)i) and Q,, := (\(/o|v(/i), an inverse 
Gram-Schmidt process yields 



Res. Ent. 



l¥o>^ = ko>^ 



(58) 
(59) 
(60) 

l¥i)^ = Q^.ko)^ + ^i-|QvPki)^, (61) 

and therefore, 

l¥r'' = ^ (ko)^ ko)^ \e,f + Q^Q^, ki)'* ko)^ ko)^ 

+ V(i-|Q^P)(i-|Q^P)ki)^ ki)'' ki)'' 



+ Q^Vi-|Q^Pki>^ki)''ko>^ 



Q^.\/l-|Q^,Pkl>^ko>^kl)^^ 

(62) 

The state is effectively a three-qubit state in 9-(2 ® ^ ® ^ ■ 
So, we can analyze the state with regard to its bipartite entan- 
glement and its residual, GHZ-like, genuine tripartite entan- 
glement Xres [40] . The relevant squared concurrences are 

C\A\B)^\Q^\\\-\(^\^) 

= |(Vo|¥i)P(l-|(^o|c^i)P) (63) 
C\A\C) = \Q^\\\~\(^\^) 

= |((^o|(^i)P(l-|(¥ol¥i)P), (64) 

(65) 



Hence, 



C2(B|C) =0, 
C\A\BC) = \-\Q^\^\%\^ 

= l-|(^o|^i)Pl(¥o|¥i)P. 



x..v = (l-|(^P)(l-|QvP) 

= (l-|(c^o|^i)P)(l-|(¥o|¥i)P), 



(66) 



(67) 



which is plotted in Fig. 10. Furthermore, the total entangle- 
"isnt C^^j^,, is given by 

QLz = C2(A|B) +C2(A|C) +C2(B|C) +T„, 



= l-|(c^o|^i)Pl(¥o|¥i) 



(68) 




Bip. Ent. 



<0O|^l> 



These results provide a basis for an interesting discussion. 
Figure 10 looks quite unspectacular. However, it contains 



Fully Sep. 



FIG. 10: (Color online) The residual entanglement of the state 
I^^ABC ^ function of the overlaps (ij)o|(|>i) and (v|/o|\|fi), which are 
assumed to be real without loss of generality. In the upper comer, 
for Trej = 1, the state is maximally GHZ-like entangled. It shows no 
bipartite entanglement. In the left and the right comers, the residual 
entanglement is zero. Nevertheless, the state is maximally entangled, 
but just bipartite entangled. Finally, in the lower corner, the state is 
fully separable. For either ((|)o|(|)i)=0 or (yo|¥i)=0' the state is al- 
ways 1-ebit entangled, while for both (())o|(|)i) 7^ and (\|/o|Yi) 7^ 
the state has less than 1-ebit entanglement. 



very useful information in combination with the other results. 
First, it can be seen that for both overlaps ((j)o|*j)i) and (\(/o|\(/i) 
being zero, the state is exactly the GHZ state, which shows 
only GHZ-like tripartite entanglement but no bipartite entan- 
glement. In the case that only one of the overlaps is zero, the 
state is always 1-ebit entangled, while it can be tuned between 
GHZ-like entanglement and common bipartite entanglement. 
If the overlap which is not zero becomes one, maximal Bell- 
state-like entanglement occurs. Finally, if none of the overlaps 
is zero the state is not 1-ebit entangled anymore and it corre- 
sponds to a mixture of bipartite entanglement and GHZ-like 
tripartite entanglement. If both overlaps are one, the state be- 
comes simply fully separable, showing no entanglement at all. 
All this is basically illustrated in Fig. 10. The upper corner 
of the graph corresponds to a GHZ state, while the left and 
the right corners represent states which are maximally bipar- 
tite entangled. Every point in the graph which is not at the 
upper left or upper right edge corresponds to a less than 1-ebit 
entangled state, while states lying on these edges are 1-ebit 
entangled. States on the lower edges of the graph are only 
bipartite entangled, while all other states show also GHZ-like 
tripartite entanglement. Finally, the lower corner corresponds 
to a fully separable state. 

As a result, just by tuning simple overlaps between the in- 
volved qumode states, it is possible to gradually switch be- 
tween these different entanglement scenarios. Such tuning of 
overlaps is a relatively easy task when the qumode states are 
simply realized by ordinary coherent states. Then only the 
amplitudes have to be adjusted. However, the experimental 
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preparation of the overall entangled tripartite state |v(/)'*^'' may 
cause difficulties. Nevertheless, this idea of tuning between 
various entanglement configurations by modification of the 
overlaps of the participating qumode states may be a scheme 
which could possibly be experimentally realized in the future. 



Non-Gaussianity of the state corresponds to non-Gaussianity 
of the characteristic function. This is equivalent to non- 
Gaussianity of the Wigner function, as (non-)Gaussian func- 
tions stay (non-)Gaussian under Fourier transformation [41]. 
Hence, consider the Wigner function of p. 



V. SUMMARY AND CONCLUSION 



W{x,p) := - I dye 



-y\p\x + y) 



In this paper, first, we proved that hybrid entanglement is 
necessarily non-Gaussian. Then, utilizing an inverse Gram- 
Schmidt process, we presented a classification scheme regard- 
ing hybrid entanglement which distinguishes between effec- 
tive DV hybrid entanglement and true hybrid entanglement. 
With aid of this characterization one can always find out 
whether a given hybrid entangled state can be analyzed by 
means of DV methods or not. To illustrate this framework, 
a few exemplary states have been discussed, and witnessing 
of true hybrid entanglement has been demonstrated exploiting 
the well-known SV inseparabiUty criteria. Finally, we exam- 
ined multipartite hybrid entanglement, especially the tripartite 
case, and briefly outlined some differences compared to the 
bipartite regime. As an example, a typical tripartite hybrid 
entangled state was investigated. 

To conclude, we want to point out that, although many ideas 
for hybrid protocols and schemes have appeared recently, no 
thorough classification of hybrid entanglement from a more 
formal point of view has been proposed so far. The present 
work is supposed to fill this gap. Nevertheless, there are 
many open questions: For example, entanglement quantifi- 
cation of true hybrid entanglement remains an unsolved prob- 
lem in general. Especially in the multipartite case, various 
questions remain. For example, what can happen when we 
consider more parties than just two or three? 
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Appendix A: PROOF OF LEMMA 1 

Lemma 1 states that any single-partite li-dimensional quan- 
tum state with finite d and d > 2 is non-Gaussian. Proof. 
Consider such a general c/-dimensional state in a pure-state 
decomposition and make use of the Fock basis of a d- 
dimensional subspace of the Fock space: 

N N 

P = E/'nIV") ' P« >OVn; ^ /7„ = 1, 

H= 1 « 

d-l d-l 

Wn) = '^nm\m) , c„,„ G CVn,m; ^ \c„,„\^ = 1. 



m=0 



: f ^ [dy^^'Py £ c„cMx^y)rk{x + y)- 

n=l kJ=0 



(A2) 



The position wavefunction of a Fock state |«) is given by 

xif„{x) ^ {x\n) ^ //„(x) ^^^^ 

where Hn{x) are the Hermite polynomials [41, 42]. Therefore, 



(.r-r)^ (.t+v)^ 



Hk{x+y)H,{x-y) 



n=\ 
d-l 



N 

E 

n=l 



(A4) 



^ V * Hkix+y)H,{x~y) 



U=0 



Now ^ / dye Py is just a Fourier transform operator 

_ 2_ 2 

with respect to y, and e^^ is a Gaussian func- 
tion. Hence, for W{x,p) being Gaussian, the sum P := 



d-l 



H^(x+y)H,{x-y) 



(Al) 



— ; , , , — also has to be Gaussian in 
both X and y. However, P is just a polynomial. Strictly speak- 
ing, it is a polynomial of finite order. Hence, P can never be 
an exponential function and therefore also not become Gaus- 
sian. There is one exception though. P may be a polynomial 
of order zero, i.e., a constant. 

So, we can assume that p is Gaussian and show that this is 
the case if and only if c„/ = V Z > 1 . is trivial, since 
for c„/ = V n , V / > 1, p = |0) (0|, which is obviously Gaus- 
sian. For we perform an induction with respect to the 
dimension d: 

Inductive Basis. For d = 2 there is only one term x^'^'^^' = 
x^ in the sum P, which comes from //^/_ i [x + y)Hii_ li^^y) ~ 
H\{x + y)H\{x — y). Its coefficient is ci^^^j /7„|c„,ip, where 
ci ^ is a real constant due to the Hermite polynomials. If 
the state is Gaussian, Y^=i Pnkn,! P = 0. Since p„>QV n and 
\cn.i I > V n, from Y.n=i Pn\cn,i P = follows c„^i = Oy n. 
Hence, c„/ = V / > 1 (/ can be only or 1 for c/ = 2) is 
proved as a necessary condition for Gaussianity for d = 2. 

Inductive Step. Assume validity of the statement for 
d — m and consider d = m + I. Then, there is only 
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one term = x^'" in the sum P, which comes 

from //(,„+ 1 ) _ 1 (x + y )//(,„+ 1 ) _ 1 (x - y ) = i/„, (x + y)H„, {x~y). 
Again with 7^ being a real constant the coefficient is 
CmL^=i Pn\c„.m\^. If the State is Gaussian, l^^^j p„|c„,„,p = 0. 
Since /?„ > V n and |c„,„,| > V n, from /?n|c„^„,p = 
follows c„ ,„ = V «. Then, from the inductive hypothesis 
it is known that c„/ = V 1 < / < — 2 = m — 1. Hence, 
c„/ = V 1 < Z < — 1 = m, which proves the validity of the 
statement for d ^m+\ and hence for all finite m. 

Concluding, a li-dimensional quantum state with finite d 
and d >2, written in the form (Al), is Gaussian if and only if 
c„i — Oy I > I, which corresponds to p = |0) (0|. Therefore, 
p is Gaussian if and only if p = |0) (0|, which is only one- 
dimensional. Hence, any t/-dimensional quantum state with 
finite d and d > 2 is non-Gaussian. ■ 



Furthermore, the binomial identity has been utilized. 

The calculation of Eq. (29) proceeds similarly. How- 
ever, instead of exploiting p^,^„„„, = 1^=0 P«' l«>£ ("I i" the 
Fock basis, we made use of the coherent state basis. The 
Glauber-Sudarshan P representation of the thermal envkon- 
ment is given by [43] 



Therefore, 



%, ,(a,a*) = — — -e 



(B8) 



Appendix B: AUXILIARY CALCULATIONS 

The derivation of Eq. (26) becomes clear, when writing the 
thermal photon noise channel with aid of an ancilla Hilbert 
space (the environment E), and joint unitary evolution on the 
overall Hilbert space. Finally, a trace operation over the an- 
cilla Hilbert space has to be performed: 



Pthermai = —, — T / d^ac |a) (a| . (B9) 



Expressing the thermal state in this form, the action of the 
thermal channel on an element |a) (p| is 



ermul) ^ 
{al-aB-agUE) 



(Bl) 



with 

(1) beam splitter unitary U^^ 

(2) environmental thermal state 

Pfhennal = ^^^=0 (i+j„'*])"+I l")^ ("I [^3], 

(3) mean thermal photon number («,/,), 

(4) and beam splitter transmissivity r| = cos^ 9. Note that 
for (nth) = 0, the photon loss channel is obtained, which is 
therefore just a limiting case of this channel. For the calcula- 
tion, the following relations have been exploited: 



^thermal 



(|a> 



^ ' d^-yt ^tr£[j7^^ (|a)^(PI 



7t («f/,; Jc 
.|y)^(y|)j7 



(BIO) 



With the relations (B2), (B4), (B5) and U'^'^D" {j)U'' 
D^(0Ty)D^(-\/r^Y) we find 



ljBE'(jBE 



-tBE' 



a 



= 1, 



:|0)^|0)^ 

:D^(a)|0)^ 



(B2) 
(B3) 
(B4) 
(B5) 
(B6) 

(B7) 



hhermali\(i)^ (PI) 



X ( Vi-iiP + \/iTyI Vi-na + \/ny> 
X l\/iTa - Vi-iiy>^ (\/iTP - Vi-iiyl ■ 

(BID 

Inserting this expression in the overall state p''*^ = (1'* (g) 
$L™«/)l¥r (¥l with = ^(|0)'^|a)V|l)^|-a)^) 
yields Eq. (29). 
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